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Abstract 

l^ We present a precise isospin analysis of the B —¥ D^D^K decays using new recent 

experimental measurements on these final states. The decays B — > D^D^K, originating 
from b — y ccs transitions, are linked by a rich set of isospin properties. The isospin 
relations that connect the decay modes are presented and a fit is performed to obtain 
the isospin amplitudes and phases. We discuss the results of the fit and present a new 
measurement of the ratio of branching fractions S(T(45) — > B + B~) and B(T(4S) — > 
B°B°). We finally discuss the implications of our findings for the measurement of the 

£Sj ! unitarity matrix parameters sin(2/3) and cos(2/3) using these decays. 
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1. Introduction 



In this Letter, we use an isospin analysis to establish relations between the different 
B — > D^D^K decays. These decays proceed via b — > ccs transitions, which are known 
to present peculiar isospin properties [l|. The possibility that a large fraction of b — > ccs 
decays hadronize as B — > D^D^K was first suggested in Ref. [2J in the context of the 
discrepancy between the measured B semi-leptonic rate and the theoretical prediction. 
This hypothesis was confirmed by many experimental results where it was found that 
B ->■ D^D^K decays account for about 4% of the B° and B + decays 0, 3, Q. These 
results provide the input for the isospin analysis and the test of the isospin relations. An 
additional motivation for an in-depth study of these channels is the possibility, originally 
discussed in Refs. p, |7|, |8(, to measure sin(2/3) and cos(2/3) using these decays. Indeed 
they proceed through the same quark current than the gold-plated mode B° —¥ J/^/K° 
and are not Cabibbo-suppressed to the difference of the B° — >• £)(*)£)(*) modes. 

This Letter, which updates and supersedes a previous investigation reported in Ref. [9| , 
presents the complete set of isospin relations for B — > D^D^K decays. They are 
compared to the measurements through a fit of the experimental data which determines 
the isospin amplitudes. There are 22 possible modes for the B — > D^*'D^*'K decays; here 
the B is either a B° or a B + , the £>M i s either a D°, D*°, D + , or D* + , the X>M is the 
charge conjugate of D^*\ and the K is either a K + or a K°. 

These decays have been the object of many experimental investigations during the past 
years. In particular, the BABAR Collaboration [5] published recently a complete set of 
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measurements of the 22 branching fractions with an excellent accuracy. They used 471 x 
10 6 BB events collected at the T(4S) resonance, corresponding to an integrated luminosity 
of 429 fb _1 . The Belle Collaboration performed a measurement of the branching fractions 
of the modes B° -»■ D*-D* + K° and B + -► D°D°K + [11] using 449 x 10 6 BB pairs. 
All these results are used in our analysis. 

With respect to the previous study |9J, the statistical and systematic precision on the 
experimental data is improved by a factor three or larger, thereby improving by the same 
amount the statistical power of the tests performed. This allows to put on a firm ground 
the conclusion that we draw from this study. 

In addition to the higher statistics, another improvement of the analysis shown in this 
Letter is the fact that the branching ratios B{T(AS) -* B+B~) and B(T(AS) ->■ B°B°), 
needed to compare the neutral to charged B meson decays measured at an e + e~ machine 
operating at the T(4S) resonance, are presently known with a good accuracy. This good 
knowledge of the ratio helps to constrain more strongly the fit performed here. 

The aim of this study is: 

• to verify the isospin relations using a new set of precise experimental results; 

• to provide some insight into the B — > D^D^K decay mechanism from the inspec- 
tion of the isospin amplitudes; 

• to discuss the implications of our findings for the measurement of sin(2/3) and cos(2/3) 
using these decays. 

It has to be noted that other authors [12J studied B — > D^D^'K decays in the context 
of color rearrangement models, and compared their predictions with the experimental 
measurements. 

2. Isospin relations for B — > D^D^K decays 

A full derivation and discussion of the isospin relations for these decays can be found 
in Ref. [9]. Here only the main results is summarized. 

The B — > D^'D^K decays proceed via a b — > ccs current through the diagrams 
of Fig. [TJ Depending on the final state, the external W-emission diagram, the internal 
W-emission diagram (which is color-suppressed), or both contribute to the transition 
amplitude. A penguin diagram, shown in Fig. [2] (left plot), can also contribute to the 
b — > ccs current. It is expected to be suppressed relatively to the tree diagrams of Fig. [I] 
and does not modify the isospin relations. 

The decays B° -* D^)°D^ K° and B + -► D^°D^°K + could also proceed through a 
different diagram, shown in Fig. [2] (right plot), which could introduce a A/ = 1 amplitude. 
However this diagram proceeds through two suppressed weak vertices b — > uW and W — > 
su and a cc pair must be extracted from the vacuum, instead of a light quark pair as 
in the Cabibbo-allowed diagrams. This amplitude is therefore suppressed by at least a 
factor A 2 , where A is the expansion parameter of the Wolfenstein parametrization. For 
these reasons we expect that AJ = holds to an excellent precision. 

As already mentioned, the isospin properties of the b — > ccs current are well known and 
follow from the fact that only isoscalar quarks are involved. Therefore this is a A/ = 
weak transition and the final state is an isospin eigenstate. 




K 



D' 



D< 




D< 



K 



{•) 



D' 



(V 



Figure 1: Left: internal VF-emission diagram for the decays B — > D^*>D^K. Right: external VF-emission 
diagram for the decays B -¥ D^D^K. 
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Figure 2: Left: QCD penguin diagram for the decays B — > D^D^*'K. Right: Cabibbo-suppressed 
diagram with A/ = 1 amplitude. 



The isospin properties translate in the following set of relations 
A ( B°^D-D+K«) = ±A l + -L Aa 
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A(B" -+ D°D°K 
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where A\ (Aq) is the amplitude to produce the system DK with an isospin quantum 
number equal to 1 (0). The Ai amplitudes in these formulae are reduced matrix elements, 
in the terms of the Wigner-Eckart theorem, of the isoscalar Hamiltonian. 
A similar set of relations holds for charged B meson decays 



A(B + -)■ D"D + K") = -^A 1 -^A Q 



A(B + ^D°D°K + ) = -—Ai + ^Ao 



i , i 

VE \/2 



(4) 
(5) 




A(B + -)• D-D + K + ) = -\ -Ax, (6) 



where the A amplitudes are the same as for the neutral B decays. Identical equations 
hold for the other set of decays, B ->■ DD*K, B -> D*DK and B -» D*D*K, with 
different amplitudes A in each case. Equivalent relations can be obtained considering 
the isospin quantum numbers of different subsystems of the final state (DD, DK). The 
DK subsystem is chosen here because in this case the transitions of Eqs. (jSJ) and (jSJ), 
proceeding only through the color-suppressed diagrams of Fig. [Q (left plot), are associated 
only to the A\ amplitude. 

The relations presented above can be cast in the form of a triangle relation between 
the amplitudes: 

- A(B° ->■ D-D°K + ) = A(B°^D-D + K°) + A(B° ^ D°D°K°) (7) 

-A{B + -> D°D + K°) = A{B + -> D°D°K + ) + A(5 + ->■ D-D + K + ), (8) 

which are depicted in Fig. |3j The two triangles for B° and i? + decays are identical 
according to the isospin relations, however experimentally it is advantageous to build the 
triangles separately with the B° and B + amplitudes. 





D" D* K*) 



A(B ° -> D" D° K*) A(B * -» D° D* K°) 



Figure 3: Isospin triangles for the B° (left) and B + (right) amplitudes. 

We finally notice that Eqs. (p^ to (JHJ) are valid not only for the total decay amplitude 
but also for each helicity amplitude separately as well as for the amplitude as a function 
of the Dalitz plot coordinates. The amplitudes and phases we extract from the fit are 
averaged over the Dalitz plot as well as over all the accessible final states (vector po- 
larizations, partial waves, ...). This remark is of particular importance since it is now 
well known that many resonances are present in the B — > D^D^K decays, such as the 
#(3770), the £> s i(2536), the X(3872), and the £> s i(2700) mesons Q. 



3. Study of the experimental results 

The branching fractions for the charged and neutral B meson decay can be written 

B(B+ -+ DUdMK) = (27r)3 T 3 2M | + (J '^W)^W) \A{B + ^D^D^K)\" (9) 
B(B°^D^D^K) = (27r)3 T 3 2M 3 o (/ 'H-WDOiHwic) \MB + -> D^ D^ K)\ 2 , (10) 



where r+ = 1.638 x 1CT 12 s and r = 1.525 x 10~ 12 s 14J are the lifetimes of the B + 
and -B mesons, M B + and M B o are the masses of the B + and B° mesons, m^) D ^) and 
m D (,) K are the invariant masses of the D^*'D^*> and D^K subsystems, and the integral 
is computed numerically over the allowed region of the three-body phase space. 

The BABAR Collaboration has recently studied the full set of B — > D^D^K decays 
and has provided precise measurements for all these modes [5j. We use also the exper- 



imental results from the Belle Collaboration [10|, LL3J which are available for the modes 
B° — y D*~ D* + K° and B + — > D°D°K + . These two modes from Belle are combined with 
the corresponding ones from BABAR assuming fully correlated systematic uncertainties. 
Table [1] presents the measurements of the B — > D^D^K final states after having com- 
bined the BABAR and Belle results. 

The BABAR and Belle data have been collected at the PEP-II and KEKB accelerators 
from the reaction e + e _ —} T(4S') — > BB. To compute the branching fractions, it has been 
assumed that B(T(AS) -+ B+B~) = B{T(AS) -* B°B°) = 0.5. However these equalities 
do not necessarily hold. In order to account for this factor, we rewrite Eqs. (J9j and ( ITU]) 
in term of the rescaled amplitudes 

A= , A (11) 

y/2B(T(AS) -> B°B°) 

The expression for B(B + — > D^D^K) is then multiplied by the additional factor 

B{T{AS) -> B+B-) 
/+/o B pr(AS) -> B°B°) ' l ' 

The experimental data are fitted simultaneously using the \ 2 method: 

( f _ /WAN 2 
X — (^cxp — ^prcd) V (B cxp — £> pre d) H 2 ' ("^J 



J +/0 

where £> exp represents the vector of the branching fraction measurements, £> pre d repre- 
sents the vector of the branching fraction predictions, and the superscript T denotes 
the transposed vector. The predictions depend on 13 parameters which are f + / and, 
for each set of decays, |Ai|, \A \, and 5 = arg(A\Ao). The matrix V is the covariance 
matrix between the 22 branching fraction measurements, which allows to take properly 
into account the systematic uncertainties that are common and correlated between each 
mode. The correlated systematic uncertainties consist of uncertainties originating from 
the signal shape, the reconstruction and the identification of particles (charged tracks, 



soft pions from D* + decays, K®, ir°, single photon, and K + identification), the branching 
fractions of the secondary decays (Z)W and Kg), and the accounting of the number of BB 
pairs produced in the experiment (see Table III of Ref. [5|). We separate each contribu- 
tion of these systematic effects in order to break down the problem into quantities which 
are completely independent or completely correlated. We sum these separate covariance 
matrices together to obtain the total covariance matrix, where the partial correlation 
structures emerge. The last term in Eq. (fT5|) constrains / + /o to the world average value 
/^ = 1.065 ± 0.026 Q. 

The results of the minimization of this x 2 are reported in Tables [1] and [2j The overall 
agreement between the measured and predicted branching fractions is fair as can be judged 
from Table (TJ Figs. H] and El and from the value \ 2 = 18-9 for 10 degrees of freedom 
(ft-dof) with a probability of 4.1%. We observe that the main source of the disagreement 
concerns the modes containing one or two D*° mesons, with a measured branching fraction 
systematically above the predicted value. This could point to a systematic shift that was 
not properly taken into account in the experimental analysis. For some B° decays which 
are not distinguishable experimentally, only the sum of the branching fraction with the 
charge conjugate final state has been measured. We present in Table [3] the fitted values 
for the individual branching fractions. 

The fit has also been conducted without the constraint on /+/o- We obtain a value 

/+/o = 1.100 ±0.056 (14) 

which is in good agreement, while less precise, with the world average. 

An alternative way of displaying the experimental results and the fit results is given 
by the isospin triangles introduced in the above. For ease of comparison, we normalize the 
triangles to the size of the basis (\A(B° -> D^~D^°K + )\ and \A(B+ -)■ D ( * )0 D w+ K°)\): 
therefore the lower side extends in each case from (0,0) to (1,0) and the shapes of the 
triangles can be directly compared. Given that we have only a measurement of the sides, 
there is a fourfold ambiguity on the vertex of the triangle. We choose consistently the 
same solution for its orientation. The seven measured triangles defined in this way are 
shown in Fig. together with the fit result. We notice that in all the cases the shape of 
the triangles presents large angles. 
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Figure 4: Results of the x 2 fit to the experimental branching fractions. The fitted branching fractions 
are shown by the stars while the points with error bars show the measured values. 
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Figure 5: Ratios of measured branching fractions over predicted branching fractions, B cxp /B prc d- The 
vertical line shows the case B cxp /B prc d = 1. 



Table 1: Branching fractions (B) for each B — > D^D^K mode. The second column shows the experi- 
mental results while the third column presents the result of the \ 2 fit- The first error on the experimental 
branching fraction is the statistical uncertainty and the second is the systematic uncertainty [51 lloL 111 | - 
The experimental results from the modes B° —> D*~ D* + K° and B + -+ D°D°K + are a combination 
between the BABAR and Belle measurements. 



B decay mode 


Bexp. (10- 4 ) 


B fit (IQ" 4 ) 


B° decays through external IV-emission amplitudes 


B° -> D-D°K + 
B° -> D-D*°K + 
B° ->■ D*-D°K + 
B° ->■ D*~D*°K + 


10.7 ±0.7 ±0.9 

34.6 ±1.8 ±3.7 

24.7 ±1.0 ±1.8 
106.0 ±3.3 ±8.6 


10.9 
31.5 
21.8 
88.4 


B° decays through external±internal iy-emission amplitudes 


B° -> D-D + K° 

B° -)■ D*-D + K° + D-D* + K° 

B° -> D*-D* + K° 


7.5 ±1.2 ±1.2 
64.1 ±3.6 ±3.9 
79.3 ±3.8 ±6.7 


10.2 
65.1 
76.7 


B° decays through internal VT-emission amplitudes 


B o ^ d°D°K° 

B° ->■ D°D*°K° + D*°D°K 

B° -»■ D*°D*°K° 


2.7 ±1.0 ±0.5 

10.8 ±3.2 ±3.6 

24 ±5.5 ±6.7 


2.0 

10.1 

10.7 


B + decays through external H^-emission amplitudes 


B + -> D°D + K° 
B + -)■ D°D* + K° 
B + _^ d*°D+K° 
B + _^ d*o D *+ K o 


15.5 ±1.7 ±1.3 
38.1 ±3.1 ±2.3 

20.6 ±3.8 ±3.0 

91.7 ±8.3 ±9.0 


12.5 
36.3 

25.1 
101.7 


B + decays through external±internal jy-emission amplitudes 


B+ -> D°D°K + 
B + ->■ D°D*°K + 
B + _^ d*°D°K+ 
B + _^ d*o B *o k + 


14.0 ±0.7 ±1.2 

63.2 ±1.9 ±4.5 

22.6 ±1.6 ±1.7 

112.3 ±3.6 ±12.6 


11.7 

55.2 
19.7 
88.3 


B + decays through internal VT-emission amplitudes 


B+ -> D-D+K+ 
B + -> D-D*+K+ 
B + _^ D*-D+K+ 
B+ -> D*-D*+K+ 


2.2 ±0.5 ±0.5 

6.3 ±0.9 ±0.6 
6.0 ±1.0 ±0.8 
13.2 ±1.3 ±1.2 


2.4 
6.0 
5.7 
12.3 



Table 2: Results of the \ 2 fit to the experimental branching fractions for the amplitudes and phases. The 
superscripts LL, L*, *L and ** refer to the B -> DDK, B ->• DD*K, B -> D*DK and B -> D*D*K 
decays respectively. The amplitude values are in units of 10 while the phases 8 are in degrees. 



Parameter 


Value 


S LL 


0.23 ±0.03 
0.59 ±0.02 

94 ±8 


\ AL*\ 

1^1 1 

1 AL*\ 

5 L * 


0.42 ± 0.04 

1.33 ±0.04 

53 ±9 


\Af\ 

A*L\ 

8* L 


0.41 ±0.04 

0.92 ±0.03 

103 ±7 


1 A ** 1 

I 71 ! 1 

1 4**1 

\ A o 1 
5** 


0.72 ± 0.05 

2.28 ±0.08 
100 ±7 


/+/0 


1.071 ±0.023 


X 2 /n do { 
Prob(x 2 ,n dof ) 


18.9/10 

4.1% 



Table 3: Fitted values of the branching fractions for the B° 
have not been measured individually. 



DD*K° and B° ->• D*DK° decays which 



B decay mode 


B fit (10- 4 ) 


B° ->■ D*-D + K° 


17.1 


B° ->■ D-D*+K° 


48.0 


B° ->• D*°D°K 


4.9 


B° ->■ D°D*°K° 


5.2 



10 









Figure 6: Isospin triangles for the B — > D^D^K amplitudes. Each plot presents the measured vertex 
of the triangle, where the basis has been normalized to unity. The dotted contour shows the one standard 
deviation region. The star shows the result of the fit. 
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4. Discussion 

4-1. Dynamical features of the amplitudes 

The amplitudes and phases extracted from the data present some distinctive fea- 
tures. First, within each set, the amplitude related to the color-suppressed decays is 
much smaller, as expected. The ratios Aq/A\ are presented in Table HI These ratios are 
very close to the naive expectation of a suppression factor N c = 3, where N c is the number 
of colors. 

Second, the central values for the relative phases 8 are in all cases large and close to 
90°. From this we can conclude that there is a firm indication for large strong phases in 
these amplitudes. This suggests the presence of non-negligible Final State Interaction for 
these decays. This is both an important indication per se and has also consequences for 
the CP violation studies that will be discussed in the next section. 

Table 4: Ratios A0/A1 from the fit to the data. The uncertainties take into account the fit correlations 
between Aq and A\. 



ratio 


value 


\A L o L \/\A^\ 


2.57 ±0.37 


\a!t\i\a?\ 


3.15 ±0.28 


\Al L \/\A\ L \ 


2.23 ±0.26 


\A* *\/\Al*\ 


3.17 ±0.21 



4-2. Implications for the measurement o/sin(2/3) and cos(2/3) 

All the B° — > D^D^K final states are in principle good candidates for the measure- 
ment of the (3 angle of the unitarity matrix p, |7J, |8|. The advantages of these modes, for 
example with respect to B° — > D^*'D^*\ are that they are Cabibbo-favored and present 
a small penguin contribution. Since both B° and B° can decay to D^*' D*-*' K° , we ex- 
pect a time-dependent CP violating asymmetry. A study of the time-dependent Dalitz 
plot allows to access the phase /3 related to the B° and B° mixing. We notice that for 
B° -> D*-D* + K°, the measured value of the branching fraction (79.3 ± 3.8 ± 6.7 x 10" 4 ) 
and the value predicted by our fit (76.7 x 10~ 4 ) are almost a factor two lower that what was 
anticipated in Ref . |8| , thereby unfortunately also reducing the comparative advantage of 
this mode with respect to B° — > D*~D* + . 

The BABAR experiment did a study of the final state B° — > D*~D* + K° in this context 
and was able to constrain cos 2/3 to be positive at the 94% confidence level (under some 
theoretical and resonant substructure assumptions, and using 230 x 10 6 BB pairs) [151 ]. 
The Belle experiment did a similar analysis on the same final state with 449 x 10 6 -B-B 
pairs and did a measurement of the CP violation parameters, although the study did not 



allow to conclude on the sign of cos 2/3 [10 



Unfortunately, up to now, no other B° — > D^*'D^K° modes have been studied in the 
context of CP violation. From the BABAR data (429 x 10 6 BB) [5], we see that the final 
state B° — > D*~D + K° + D~D* + K° is observed with a significance of 13<r, where a is the 
standard deviation, which shows that a CP-violation analysis would be possible. 
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For B° —7- D~D + K°, a value of 7.5±1.2±1.2xl0~ 4 is reported (with a 5a significance). 
In this case too, the estimated value of Ref. [7J (90 x 10~ 4 ) is a factor 12 above the 
measurement. However, we stress that this channel is a good candidate for CP- violation 
studies because of the nature of the final state with three pseudoscalar particles. This 
will facilitate the angular analysis to determine the helicity amplitudes. 

Finally we notice that the B° ->■ D*~D + K° and B° ->■ D~D* + K° decay modes lead 
to final states accessible to both B° and B°. They can therefore be analyzed in the same 



way as described in Ref. 16[. The strong phases play an important role for this analysis 
as the time-dependent CP-asymmetry amplitudes are proportional to sin(2/3±5'), where 
5' is the strong phase difference between A(B° -> D~D* + K°) and A(B° ->■ D-D* + K°). 
The possibly large values of the strong phases noticed in the above need to be taken into 
account for any estimate of the sensitivities of this analysis. 

5. Conclusion 

We have presented an isospin analysis of the B —$■ D^*'D^*'K decays, based on recent 
and precise measurements of these final states. A fit was performed using the isospin 
relations between the different final states. We find a good agreement between the ex- 
perimental values and the fitted values. The isospin amplitudes exhibit several peculiar 
features like the presence of color-suppression and large relative phases. We find a value 
of i/ r( - 49 V' gogh equal to 1.100 ± 0.056, in agreement with other determinations of this 
quantity. We have discussed the features of our result and showed the implications for 
CP-violation measurements using these decays. 
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